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o ■ Abstract 

■ 

p • In the framework of percolation of color sources, the transverse momentum distribution 

D I in heavy ion and p + p collisions at all centralities and energies are shown to follow 

; 

>• ■ a universal behaviour. The width of the distribution depends on the fluctuations of 

><■ 

^ ■ the number of color sources per cluster. At low densities, there are only independent 

! 

single color sources, no fluctuations occurs and the distribution is described by a single 
exponential. At very high densities, only one cluster of many color sources appears and 
therefore there are not fluctuations either and the hardness of the distribution is sup- 
pressed. The Cronin effect in this framework is due to a maximum of the fluctuations 
which decreases as the density increases. We obtain a good agreement with exper- 
imental data including the low px behaviour and the spectra for different particles. 
We show that the transverse momentum and multiplicity distributions are related to 
each other in a deflned way. This point is satisfled by the experimental data on p + p 
collisions at different energies. 



I. INTRODUCTION 



Many data have already been collected and analized during past few years at the 
Relativistic Heavy Ion Collider (RHIC) in order to obtain a complete understanding 
of the dense QCD matter which is created in high energy heavy ion collisions. 

These data show that the inclusive high px hadron production in Au + Au central 
collisions is strongly suppressed [1]- [3] compared to the scaling with the number of 
binary nucleon-nucleon collision, Ncoiu expected on the basis of the factorization theo- 
rem for hard processes in perturbative QCD(pQCD) [4]. The suppression is larger at 
forward than at mid rapidities [5]. No suppression was found in d + Au collisions [6]- 
[7] at mid rapidity. Furthermore, the proton and antiproton yields become similar to 
the pion one at ~ 1 — 2 GeV/ c [6] . The ratio between the yields from central and 
peripheral Au + Au collisions is larger for protons than for pions at pr ~ 2 — 5 GeV/ c 
[6]. On the other hand, at low pt, the yield for pions is larger than the one for kaons, 
and both are larger than the one for protons [8]- [9]. At very low pt, the spectra for 
all species show a characteristic behaviour [9] . 

The data also show disappearance of back to back jet-like hadron correlations in 
Au + Au collisions, contrary to what is observed in d + Au and p + p collisions [10]- [11] 
and a peculiar behaviour of the fluctuations in transverse momentum [12]- [13], with a 
maximum at a certain centrality given by a number of participants around Npart ~ 150. 

All these data have originated a lot of discussion and different explanations, aimed 
to discriminate which of these effects are caused by initial state interactions and which 
ones are a consequence of final state interactions. 

Here we arc going to concentrate on the explanations based on initial state inter- 
actions. Different mechanisms of initial state interactions have been proposed, like 
saturation in high density QCD through the Colour Glass Condensate (CGC) [14]- [16] 
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or clustering of strings [17] in the string models. Another possibility is the shadowing 
through pomeron interaction in the Dual Parton Model [18]. All these mechanisms 
have in common the modification of a multiple scattering pattern -in the target rest 
frame- or gluon interaction -in a fast moving frame-. 

In the framework of the Colour Glass Condensate, the number of gluons in the 
hadron wave function can reach saturation when their momenta become smaller than a 
saturation scale Qs{x) ~ x~^^A^/^. Clustering and percolation of color sources, strings 
or partons, can also be seen as an initial state interaction phenomenon. In the case of 
strings, these are formed by color fields stretched between partons of the projectile and 
the target, located at the ends of the strings. These partons are embedded in the wave 
function of the projectile and the target [19]. In this case, above a critical density, a 
cluster of color sources -strings- is formed through the whole collision area. 

In both approaches, if the number of partons exceeds a critical quantity, they will 
begin to overlap in the transverse plane, and interact with each other, which prevents 
further growth of parton densities. As the number of strings is mainly determined by 
the number of inelastic parton-parton collisions, the density of strings is not a property 
of the isolated projectile. Nevertheless, since the density of strings is connected with 
the density of partons and since the interaction and percolation of strings takes place 
before the formation of final secondaries -before fragmentation-, we call our effect 
initial. 

From this point of view, the suppression of pr and the reduction of multiplicities 
have the same origin [20]- [21], namely, the overlapping of color sources. In both 
approaches, the transverse momentum distributions satisfy a scaling law which is in 
agreement with the experimental data. In the CGC, the suppression was predicted 
to be stronger outside mid rapidity and this point was later confirmed [5]. Below, we 
will show that this is also the case in the framework of percolation of strings. A more 
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quantitative comparison of the results of both approaches shows also a remarkable 
agreement [22]. The fact that the results of the string clustering approach, which is a 
soft, QCD inspired but model dependent description, coincide with the results of the 
CGC picture, which is a theory deduced from perturbative QCD, induce us to think 
about the possibility of a smooth transition between soft and hard regimes and about 
a perturbative confinement [23]. 

All these features, together with a correct description of the transverse momentum 
fluctuations and of the low pr data, point out the existence of initial state interactions. 
On the other hand, the disappearance of back to back jet-like hadron correlations in 
Au + Au collisions, contrary to what happens in d + Au collisions, seems to point 
out a final state interaction description such as jet-quenching [24] or interactions with 
partons and comovers [25]. 

The Cronin effect is crucial in order to discriminate which kind of interaction is 
working [26] . In this paper we will show that this effect is a low energy one, and it is 
going to disappear at higher energies and/or densities. Similar results are obtained in 
the framework of the CGC [27]. 

Before starting with the description of our model, let us remember some general 
features concerning the clustering. On general grounds, in the clustering approach, one 
has to distinguish between two density regimes, high and low. Let us give a very simple 
example: the problem of throwing Ns coins into M boxes (one can read: production 
of Ns partons or strings in the interaction area ttR"^). The distribution P{N) of N 
coins in a box (one can read: N partons or strings in a cluster) can be studied and, in 
particular, the inverse normalized fluctuation k 

<N^>-<N>^ ^ ^ 

can be analytically calculated. At small density, 77 = ^ << 1, the coins are isolated 
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and A; — > oo. At large density, rj » 1, the coins are equally distributed in the boxes 
and, again A; — > oo. At some intermediate value rjmin, there is a minimum. In the low 
density regime rj < rjmin, k decreases as 77 increases, while in the high density regime 
V > Vmin, k increases with rj. 

This behaviour for the fluctuations in the number of partons or strings in a cluster 
will be crucial in order to explain most of the data. In this paper we will explore these 
points in the framework of string percolation, using the strings as our basic objects. 
The use of strings could be considered as very model dependent. Nevertheless, the 
string structure can be derived directly from QCD under certain approximations [28]- 
[30]. In fact, a string profile can be derived from the study of correlators of gluons 
fields, profile which is in agreement with the results obtained independently in lattice 
QCD [30]. 

In order to apply these ideas, we will use the following ingredients: We need to know 
the number of strings Ns- Up to RHIC energies, Ns in the central rapidity region is 
approximately twice the number of collisions, A^coii- However Ns can be larger at RHIC 
and LHC energies. We compute Ns using the Quark-Gluon String Model, equivalent 
to the Dual Parton Model [31]- [32]. Most of the reasonable string models [33]- [39] of 
heavy ions collisions obtained similar results for A^^. This fact gives us confidence in our 
values. Notice that sometimes, even in experimental analyses, Ncoii is obtained from 
the Glauber model without taking into account the energy-momentum conservation. 
At high energy this conservation reduces Ncou- In our computation of Ns we use 
a Monte-Carlo code [39], based on the Quark-Gluon String Model, which takes into 
account energy-momentum conservation. 

In order to obtain analytical formulae which can give an insight of the physical 
grounds we use soft {exp{—ap^)) transverse momentum distribution [40] and Poisson- 
like multiplicity distribution for the fragmentation of one string. We are aware that 
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a string can also produce hard particles as in the Lund string model [41], and our 
simplification must be seen as a first approximation which clearly would fail at very 
high pt- In spite of this, we will show how the clustering of strings gives rise to a uni- 
versal behaviour of both the transverse momentum distribution and the multiplicity 
distribution, which are related to each other through a gamma distribution that repre- 
sents the cluster weight function. By the overlapping of soft strings, the soft spectrum 
gets a hard-like contribution. This behaviour must be seen as complementary to what 
happens in parton saturation in the CGC, where the clustering of gluons results in a 
softer spectrum. At RHIC and LHC energies, there is a large range of pt where both 
descriptions should coincide. Of course, we do not claim to give a description for the 
whole Pt range at fixed energy. Our explanation applies from low to intermediate pt 
range. This range increases with the energy. At very high p^, a perturbative QCD 
description will be necessary. 

The plan of the paper is as follows: first we describe our approach and we derive 
the transverse momentum and multiplicity distributions, in order to compare them 
with experimental data in the next chapter. The comparison includes p-\- p data and 
predictions for LHC energies. After, we discuss the effect of the clustering on the 
disappearance of back to back jet-like hadron correlations. We will finish with some 
conclusions. 

II. PERCOLATION OF STRINGS, TRANSVERSE MOMENTUM AND 
MULTIPLICITY DISTRIBUTIONS 

Multiparticle production is currently described in terms of color strings stretched 
between the partons of the projectile and the target. These strings decay into new ones 
by sea q — q production, and subsequently hadronize to produce the observed hadrons. 
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The color in these strings is confined to a small area in the transverse space, Trr^, 
with tq ~ 0.2 — 0.25 fm. This value is obtained in the vacuum correlator method and 
corresponds to the correlation length of the QCD vacuum. This value is in accordance 
with lattice results. 

With increasing energy and/or atomic number of the colliding particles, the number 
of exchanged strings grows, and they start to overlap, forming clusters, very much like 
disks in the continuum two-dimensional percolation theory. At a certain critical density 
rjc — 1.18 — 1.5 (depending on the type of employed profile functions -homogeneous 
or Wood-Saxon-) a macroscopical cluster appears, which marks the percolation phase 
transition. For nuclear collisions, this density corresponds to 77 = where Ns is 

the total number of strings created in the collision, each one of an area Si — irr^ and 
Sa corresponds to the nuclear overlapping area so it depends on the impact parameter 
b. For very central collisions b — and Sa — ttRa- 

The percolation theory governs the geometrical pattern of the string clustering. Its 
observable implications, however, require the introduction of some dynamics in order 
to describe the behaviour of the cluster formed by several overlapping strings. We 
assume that a cluster of n strings behaves as a single string with an energy-momentum 
that corresponds to the sum of the energy-momenta of the overlapping strings [17], 
[20], and with a higher color field, corresponding to the vectorial sum of the color 
charges of each individual Qi string. The resulting color field covers the area Sn of 
the cluster. As Qn= (J2i QiY, and the individual string colors may be arbitrarily 
oriented, the average is zero, so Q^— nQ\. depends also on the area 5"! of 

each individual string that comes into the cluster, as well as on the total area of the 
cluster Sr, 




(2) 
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We take Si constant and equal to a disc of radius tq. Sn corresponds to the total 
area occupied by n discs, which of course can be different for different configurations 
even if the clusters have the same number of strings. One could make reasonable 
alternative assumptions about the interaction among the strings, as it was studied 
previously [20] , but the comparison with the data on the dependence of the strength of 
the two-body [42] and three-body Bose- Einstein correlations of the multiphcities [43] , 
clearly favours (2). 

Notice that if the strings are just touching each other, Sn = nSi and Qn = nQi, so 
the strings act independently to each other. On the contrary, if they fully overlap Sn — 
Si and Qn = y/nQi, so we obtain a reduction of the color charge. Knowing the color 
charge Qn, one can compute the multiplicity and the mean transverse momentum 
squared < of the particles produced by a cluster, which are proportional to the 

color charge and color field [20]- [21], respectively 



nSn 2 nSi 2 /o\ 

/^n = y^/^l , <PT>n^d-g- <Pt>1 (3) 

where /ii and < >i are the mean multiplicity and mean of particles produced 
by a single string. We observe 



^In < Pt >n= Wl < Pt >l , , ^2 ^ =^ ^^2 ^ " 

^ Pt >~>l ^ Pt ^ I 



The first relation denotes that the product is an extensive quantity, while the second 
one indicates that each cluster satisfies a scaling law that it is nothing but the Gauss 
Theorem. Prom the Schwinger formula, one obtains /ii = Si < p\ >i. Prom the 
experimental data we can fix /xi and < p^ >i, obtaining tq ~ 0.2 — 0.25 fm in agreement 
with the above mentioned QCD result. 

Moreover, in the limit of high density 77, one obtains 

Si_ ^ rj ^ 1 , . 

^""Sn"^' l-expi-T))- F{r)y ^ ^ 
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and the eqs. (3) transform into the analytical ones [20] 



A* = N strings F{vi) Hi 



<pI> 



1 



<Pt >i 



(6) 



F{r]) 



where and < > are the total multiplicity and mean transverse momentum and 
Ngtrings IS the total number of created strings in the considered rapidity range. 

In the mid rapidity region, the number of strings NgtHngs is proportional to the 
number of ^4 + ^4 collisions, N^ou ~ N'^^ [44]- [45], being A^^ the number of wounded 
nucleons of one nucleus. In this case, the density of strings becomes rj — N strings^ ~ 
NjI . At high densities, one should consider the percolation limit -all the strings overlap 
into a single cluster that occupies the whole nuclear overlap are-. In this case, from eq. 
(3) one obtains 



In other words, the multiplicity per participant becomes independent of Na-, i.e. satu- 
rates. 

Outside mid rapidity, Ngtrings is proportional to the number of participants Na 
instead of to the number of collisions Njl . Therefore, there is an additional suppres- 
sion factor N]1^ compared to central rapidity. This fact is at the origin of the larger 
suppression at y = 3 of the observed pt distributions in Au-\- Au and d-\- Au collisions. 

We use eq. (3) to compute the multiphcities, using a Monte-Carlo code based on 
the Quark-Gluon String Model to generate the strings [39]. Each string is produced 
at an identified impact parameter. From this, knowing the transverse area of each 
string, we identify all the clusters formed in each collisions and subsequently compute 
for each of them its multiplicity in units of /xi. The value of //i was fixed by the 
comparison of our results with WA 98 data for Ph -\- Pb central coUisions. Our results 
are in agreement with SPS and RHIC multiplicity data [21]. Using the first of the eqs. 
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(6) we obtain very similar results [45]. We observe a weak dependence on Na of the 
rapidity density per participant at high centrality, i.e. saturation. On the other hand, 
in the fragmentation region, we expect the particle density per participant nucleon to 
be equal or even less than the nucleon-nucleon rapidity density [45]. Both features, 
saturation and fragmentation scaling, are in agreement with experimental data. 

Until now we have presented our results concerning mean values -mean multiplicity 
and mean pt-- In order to get distributions, we will develop the following strategy: 
we will introduce the multiplicity and the transverse momentum distribution for the 
fragmentation of one string. We will weight them with the cluster function, and then 
we will introduce the effect of the overlapping through F{r]). We will check that we 
obtain the relation (6) for the mean values. 

The multiplicity distribution in heavy ion collisions can be expressed [46] as a 
superposition of Poisson distributions with different mean multiplicities. 



The Poisson distribution P{N,n) — ^ N —< n >, represents the cluster frag- 

mentation, while the weight factor W{N) reflects the mean multiplicity distribution 
of the clusters, related to the cluster size distribution and to the number of strings 
per cluster. This quantity has contributions due to both the nuclear structure and the 
parton distribution inside the nucleon. 

Concerning the transverse momentum distribution, one needs the distribution 
/{xjPt) for each string or cluster, and the mean squared transverse momentum dis- 
tribution of the clusters, W{x), which is also related to the cluster size distribu- 
tion through the cluster tension. For f{x,pT) we assume the Schwinger formula, 
f{x,pT) = exp{—p'^x), used also for the fragmentation of a Lund string. In this 
formula x is related to the string tension, or equivalently to the mean transverse mass 




(7) 
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of the string. Assuming that a cluster behaves similarly to a single string but with dif- 
ferent string tension, that depends on the number of strings that come into the cluster, 
we can write for the total pr distribution 

f{pt) = / W{x) f{x,pT) . (8) 

The weight function W{x) obeys the gamma distribution 

W{x)^^^{^xf-' exp(-7x) . (9) 

The reason to chose a gamma distribution is the following: in peripheral heavy ion 
collisions, the density of strings is small and therefore there is no overlapping. The 
cluster size distribution in this case is peaked at low values. As the centrality increases, 
the density of strings also increases, so there is more and more overlapping among the 
strings. The cluster size distribution is strongly modified, according to the picture 
shown in Fig. 1, where we have plotted three cluster distributions that correspond to 
three increasing centralities of the collision. Each curve in Fig. 1 can be compared to 
a gamma distribution, with different k values. Moreover, the increase of centrality can 
be seen as a transformation of the cluster size distribution of the type 

^ > <N> <Nk> ^ ' 

This kind of transformations were studied long ago by Jona-Lasinio in connection to 
the renormalization group in probabilistic theory [47]. Actually, an increase of the 
centrality implies a transformation from clusters with very few strings to another set 
of cluster with a higher number of strings, as can be seen in Fig. 1, and a renormaliza- 
tion of the main variables of the clusters, i.e. mean transverse momentum and mean 
multiplicity, induced by the higher color of the new clusters. These transformations 
have been used to study the probability distribution associated to rare events [48]- [49] . 
Introducing (9) into (7) and (8) we obtain 
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r(n + l)r(A;) (1 + 7')"+*^ Jo n\ r{k) 

and 

1^=/ (ixexp(-p^x)— — (7x) -^exp(-7x) . (12) 

(1 _|_ Pr^k Jo i (/Cj 

The distribution obtained in (11) is the well known negative binomial distribution, 
whose mean value and dispersion verify 

k < > - < N 1 <n^>-<n>^ 1 1 

<n>=<N>^- — = = - + .13 

Y < N >^ k <n>^ k <n> 

In the distribution (12) the corresponding values are: 

k < > - <x >^ 1 ,^ ^, 

< X >= - , = - . 14) 

7 <x >^ k 

The parameters 7 and 7' are different, since < N >^< x >, while k is the same in 
both equations. 

The equations (11) and (12) can be seen as a superposition of chaotic sources - 
clusters- where | fixes the transverse momentum fluctuations. At small density, 77 << 
1, the strings are isolated and A; — > 00. When the density increases, there will be 
some overlapping of strings forming clusters, increasing the denominator of eq. (1) and 
therefore decreasing k. The minimum of k will be reached when the fluctuations in the 
number of strings per cluster reach its maximum. Above this point, increasing 77, these 
fluctuations decrease and k increases. 

The distribution W{x) satisfies the Koba-Nielsen-Olsen scahng, i. e. < x > W{x) 
depends only on ^f^*- This property [50] stems from the type of the transformations 
(10). 



*However, it can depend on the energy through k. 
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In order to take into account the effect of overlapping of strings on the multi- 
phcity and the mean transverse momentum, we need to included our factor F{7]) 
in the corresponding fragmentation functions, so the Poisson distribution transforms 

into P{N,n) = ^"^^(^^iv)) g^j^^ ^j-^g Schwinger formula transforms into f(x,mT) = 

exp(— m|.xF(77)). Indeed, what is happening is that the invariance of the weight func- 
tion under the transformation x ^ Xx and 7 7/ A (or equivalently N XN 
and 7' — > ^' /X), where A = F{r]), leads to the changes < p\ >— >< p\ > /X and 
<n>^A<n>in the transverse mass and multiplicity distributions respectively. 

The multiplicity distribution becomes then the universal function 
V{n + k) {i/F{r])f V{n + k) {k/ {< n >i F{ri))f 



T{n + l)r(A;) (1 + {i/F{r]))Y+^ T{n + l)V{k) (1 + ((&/(< n >i F(r/))))"+'= 
and the transverse momentum distribution behaves as 



(15) 



dN k-l 1 k-l 1 



The above equation has been normalized to 



dN 
dp\dy 



fiPT,y) = -rT:r.- (17) 



In the new distributions, < p"^ >=< p^ >oid /F{rj) and < n >= F{rj) < n >oid, 
compared to eqs. (11) and (12). This agrees with our result from eqs. (6). 
We have, in particular 

<n>= Fir^)N, < n >,= F{r^)r^{^f <n>„ < pi >= ^ ^ ^ ^ f^^' (18) 
and the p^ dispersion 



The eqs. (15) and (16) summarize our main results. The multiplicity and transverse 
momentum distributions for any type of collision and degree of centrality are universal 
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functions which depend only on one parameter, < ni > and < >i respectively. 
These parameters are related to 7' and 7 through eqs. (13) and (14). Because of this, 
our parameters 7 and 7' in the weight function are different. The additional parameter 
k depends on rj in the way that has been pointed out before. It is related to the 
fluctuations in the number of strings per cluster by eq. (1) and plays an important 
role in the behaviour of the dependence of the transverse momentum fluctuations on 
the number of participants [51]- [52]. Notice that both distributions are obtained 
from the same kernel, the gamma distribution, with the same parameter k. This fact 
implies that they are related to each other. In particular, the suppression of high pt 
production in (16) is controlled by k which also provides us the width of the multiplicity 
distribution at high multiplicities. Due to this, we predict that in p + p collisions at 
LHC energies, where the string densities will be very high, the k value will be higher 
and in consequence both distributions, the multiplicity and transverse momentum one, 
will be narrower. 

The relation between both distributions can be checked directly with the existing 
experimental data from p + p collisions. In fact, fixing the /c-values from fits to the 
multiplicity distributions, we have fitted the transverse momentum distributions, with 
the power law A{1 + bp^)"'^ obtaining an overall agreement with experimental data at 

— 23, 200, 630 and 1800 GeV. The detailed comparison is done in reference [53]. 

The values < n >i and < p"^ >i stand for the average over all particle. If we are 
interested on the study of a particular particle specie i, we shall write < n >ii and 
< Pt >ii- 

At f] —>■ 00, k ^ 00 and the distribution (16) becomes exp{—F{r])p'^/ < p'j. >ij) 
very similar to the behaviour at — > 0, exp(— < p^ >ii). 
From (16) one can calculate 
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As 0, this reduces to 

-2F{rj)pl/ <pl>ii, (21) 

while for large pr, it becomes —2k. 

We use a gamma distribution for the clusters, which is reasonable due to the above 
mentioned arguments. It may be possible to use other distributions which have the 
same properties concerning the dependence of their mean value and dispersion on 
the centrality. The resulting multiplicity and transverse momentum distributions in 
that case would be very similar to (15) and (16). Nevertheless, we have additional 
confidence in our assumption, since the resulting multiplicity and transverse momentum 
distributions have a reasonable agreement with experimental data. In fact, the negative 
binomial distribution is in good agreement with the multiplicity distribution for p + p 
collisions in a broad range of energies, and it can also describe the h + A and A + A 
multiplicity distributions. 

The power-like behaviour (p^)"'^ found for the transverse momentum distribution, 
with an exponent related to some intrinsic fluctuations, is common to many apparently 
different systems, as sociological, biological or informatics ones [54]- [56]. Distributions 
that describe the growth of the wealth of people living in stable economical systems, 
the distribution of the citations of the scientific works, or other complex networks where 
the probability P(m) of having a given node with m links is described by the free scale 
power law P{m) ~ (m)"*^ with k close to 3 obey the same behaviour. Also, it has been 
shown [56]- [57] that maximization of the non extensive information Tsallis entropy 
leads to the distribution (16). 

This universal behaviour indicates the importance of the common features present 
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in those phenomena, namely, the cluster structure and the fluctuations in the number 
of objects per cluster. 



III. DENSITY DEPENDENCE OF THE MEAN TRANSVERSE 

MOMENTUM AND K 

A more specific test of these ideas can be made in a straightforward manner from 
eqs. (18) and (19) and eq. (4): there is a universal relation between 



<Pj> 
<P^>i 



and 



1 dn 

A 



<Pt> 
^ < Pt >i 



k 



1 dn 



(22) 



Mi<n>iJ ]l k-2 F{r))^\^ N^f dy ' 

A similar relation was suggested in the framework of the CGC model [58], and a 
quantitative discussion was presented in [59]. One should notice that (22) is fairly 
independent of the kind of particle i (tt, k, p), independent of energy (from 200 GeV 
to 1800 GeV), and independent of the interacting nuclei {Au + Au at RHIC, p + p at 
Tevatron) . 

Prom comparison with data, it can be checked that k{r]) is a function with the 
properties enunciated in the introduction. 



k{ri) 



(23) 



with 



(24) 



The function k{r]) has the qualitative behaviour that we expect, namely, it has a 
minimum kmin — 3 at rjmin — 2 and goes to infinite as r; — > and rj ^ oo. This 
means that in such limits only one kind of cluster is produced. However, the minimum 
value of k is reached at a slightly high value of rj. There are reasons to think about 
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a modification of (24), keeping the right hmits for A; at 77 —> and 77 — > 00, but 
introducing a shift in rj. In fact, F{rj)^ = 1 — e"'' is the average fraction of the total 
area occupied by clusters at density 77 when the area is homogeneous. Actually, in a 
heavy ion collisions, the surface is not homogeneous when nuclei profile functions of 
the Wood-Saxons type are taken into account [60] . This increases the average fraction 
of occupied area, resulting in a faster decrease of the exponential. Due to this, we use 
the function given by eqs. (23) and (24) for k but introducing the shift 77 — > 677, in 
such a way that the minimum value of k is reached at T/^m = 0.6 (6 = 3.3). In Fig. 
2 we present the functions before and after the change. Notice that in both cases the 
minimum value is ~ 3, where the variance of the p\ distribution diverges. 

The main dependence of k on the energy comes from the dependence on 77. In 
A-\-A coUisions, as the energy increases, the density of strings increases and k increases. 
However, there is an additional dependence on the energy, working on the opposite way, 
which up to now has not been taken into account. It is well known that even in hadron- 
hadron collisions at ^/s ~ 30 GeV , there is enough energy to produce hard scatterings. 
This modifies the single exponential in the transverse momentum distribution. As 
far as the number of strings is very low there is no possibility of overlapping. In 
perturbative QCD, the study of the production of n gluons leads to a distribution 
whose width is controlled by /c/i, being /i = 1/(1 — 7), 7 = ~ + ^)^- 
the energy increases, the QCD coupling constant c^s — > and // — > 1, in such a way 
that asymptotically only the parameter k remains. As kii decreases with the energy 
the distribution becomes broader. 

To take into account this energy dependence we will use an effective k given by 

k,ff ^k[a{^x-l) + l] (25) 

with // as defined above. We use as{Mz) — 0.122, so as{\/s — 19.4GeV) ^ 0.165, 
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a^{^ = 200GeV) ^ 0.108 and as{^/s = 1800GeV) ^ 0.081. The above expression 
satisfies kgff = /c for s — >• cxd. The values for as are lower that the ones used in other 
approaches for hadronic interactions, especially at SPS energies -at low energy, thresh- 
old effects can be not negligible-. In other to take into account these uncertainties, 
we have introduce the phenomenological expression (25). In this way, the dependence 
of k on CKg is modified through the factor a. The value of a, a = 0.3, was determined 
by the comparison of the pr distributions for Pb + Pb collisions at SPS energies and 
Au + Au collisions at ^/s — 200 GeV with the experimental data. Prom now on, in 
order to compare with experimental data, we will use k^ff from eq. (25). Through 
this equation, k -solid line in Fig. 2- is multiplied by a factor that depends on the 
energy. This factor is around 1.28 at SPS energies, 1.19 at RHIC energies and 1.13 
at LHC energies (5500 GeV). Note that the dependece of A; on 77 and on the energy is 
model dependent, so it may suffer changes, specially at low energy where the threshold 
effects can be not negligible. The fact that we want to stress here is the existence of 
a minimum on the behaviour of k versus 77, related to a maximun of fluctuations, and 
the limits of k at low and high density, where /c — > 00. 

IV. COMPARISON WITH THE EXPERIMENTAL DATA 

As we said before, the number of strings is computed using a Monte-Carlo code [39] 
based on the quark gluon string model. Knowing Ns and Tq — 0.25 fm, we can compute 
rj for each type of collision. Prom eqs. (23)-(25) we obtain k. Pinally, the values of 
< Pt >ip— 0.30, < >ik— 0.14, < p"^ >i7r= 0.06 are normalized to reproduce the 
dependence of < pf. > on centrality for protons, kaons and pions. Note that the 
different < p^ >is correspond to different 7s in the original gamma distribution 

With these values, we have entirely determined eq. (16) for all type of coUision 
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energies and rapidities. 

At low pt-i the behaviour given by eq. (21) is clearly confirmed by PHOBOS data, 
namely the derivative of log/(pT) with logpT vanishes as pt — > 0. As < >ip > 
< Pt > < Pt the absolute value is larger for pions than for kaons and than 
for protons. At higher pr, the distributions become similar. This is shown in Fig. 3, 
where we present our results for Au + Au central coUisions at yfs — 200 GeV, r] = 2.5, 
keff — 4, together with the PHOBOS data at low pt and PHENIX data at higher pt- 
Notice that in the figure there are data from two different collaborations with different 
normalizations. Probably, this is at the origin of the minor differences between data 
and our results. (We do not fit the PHOBOS data, but just apply equation (16) fixing 
the normalization ^ to reproduce the point at pt — 2 GeV) . 
Now, let us discuss the interplay between low and high pt- 
One defines the ratio Rcp{pt) between central and peripheral collisions as 

iicPyPT) - -J] TTTTTT l^Oj 

where the distribution in the numerator corresponds to higher densities, rj' > rj. The 

dn 
dy 



division by Ncou essentially eliminates Ns from ^ (this is true only at mid rapidity) 



and from (16) and (6) we obtain 



'2 ^ ^ 



') Pt 



<pf >ii 

In the Pt ^ hmit, taking into account that | < ^ < 1 and that F{ri') < F{ri), we 
obtain 

^'CP(O) - < 1 , (28) 

approximately independent of k and k'. As r]'/r] increases, the ratio Rcp decreases, in 
agreement with experimental data. 
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As pt increases, we have 



Rcp{Pt) 



1 + F{r])pl/ < >u 
1 + F{7i')p't/ <PT>li ' 



(29) 



and Rcp increases with pt (again, F{rj) > F{r)')). 



At large pt, 



Rcp{pt) 




F{r]) k' 2k-k' 



(30) 



which means that if we are in the low density (low energy, low Na) branch of the k{r)) 
curve, see Fig. 2, k > k' and Rcp{pt) > 1) (Cronin effect). As rj'/r] increases the ratio 
Rcp{pt oo) increases (it must have a limit due to phase space limitations). 

As we increase the energy of the nucleus-nucleus collision, the energy density in- 
creases and one has to observe the high density branch of the k{'q) curve. There, for 
T]' > Tj, k' > k and suppression on pt occurs. The Cronin effect will disappear at high 
energies and/or densities. 

In the forward rapidity region, the division by Ncoii in (26) does not cancel Ns from 
^, since in this rapidity region Ns is proportional to Na instead of Ncou- Therefore, 

AT' /TV' 

an additional factor jv^/jv°°il appears now in Rcp{pt)- As N'^^^ — N'^ is much larger 
than Ncoii — Na, Rcpipr, y = 3) < Rcpipr, U — 0), thus a further suppression occurs, 
in agreement with experimental data [61], [62]. In order to quantify this suppression, 
let us consider central and peripheral d + Au collisions at mid rapidity and at y = 
3.2. Prom the BRAHMS Collaboration data [61] at =^ 2 - 3 GeV/c, Rcp{pT,y = 
3.2) /Rcp{pt, 2/ = 0) ~ 0.45/1.25 = 0.36 for - 20% and 60 - 80% degree of centrality. 
For these centralities, the quoted values by the collaboration are N'^{d) — 1.96, N'^^^ = 
13.6, NA{d) = 1.39 and N^u = 3.3. The resulting ratio (1.96/13.6)/(1.39/3.3) = 0.35 
is in perfect agreement with data. 

The ratio between Rcpipr) for two different particles, for instance p and tt, becomes, 
at large px, 
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( 



< Pt >i^ 



) 



k'-k 



(31) 



For Au + Au central and peripheral collisions at = 200 GeV, we have k' — 4 
and k — 3.6 respectively. Therefore, the ratio is close to 2, in good agreement with 
experimental data. 

In Fig. 4, we compare our results for 7r° production in central and peripheral 
Au + Au and p + p collisions at ^/s = 200 GeV, together with the experimental data. 
In Fig. 5, the nuclear modification factor for central and peripheral Au + Au collisions 
at mid rapidity are shown. In Fig. 6, we show the comparison of our results for 
d + Au coUisions at ^/s — 200 GeV and at midrapidity, and in Fig. 7 we show the 
corresponding nuclear modification factor. In Fig. 8 we compare our results with 
experimental data [63]^ for central and peripheral Pb + Pb coUisions at SPS energies. 
We observe a reasonable agreement in all cases. Notice that the trend of data can 
be understood qualitatively from our eq. (16) and the ratios Rcp{pt) or the nuclear 
modification factor depend on the difference k'^r]') — k{r]) (at large pt) and on the ratio 
(at low Pt), and not on the absolute values of k', k, F{r)') and F{r)), in such a way 
that the uncertainties in the computation of 77' and 77 due to the values of Tq and Ns 
are essentially canceled. All the explanation stems from the facts that ^(77') < F{r]) 
for Tj' > rj and k' > k for large densities and k' < k for small densities. Finally, in Fig. 
9, we show our prediction for LHC energies. 

'''We are aware that there is a big uncertainty concerning SPS data (see Ref. [64] for more 
details). Nevertheless, we use these data since they have been obtained by the same experi- 
mental collaboration as the ones we have used to fix the value of /xi in our model. With this 
value we can reproduce also multiplicities at SPS and RHIC energy. These considerations let 
us to have a consistent description for pT and multiplicity distributions. 
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V. BACK JET-LIKE HADRON CORRELATIONS 



One of the most interesting data from RHIC is the disappearance of back to back 
jet-hke hadrons in Au + Au coUisions, contrary to what happens in o? + Au coUisions. 
We think that the suppression of produced jets in central A-\-A coUision is indeed due 
to final state interaction -jet-quenching or interactions with partons and comovers-. 
Here, for completeness, we would like to evaluate the interaction of a quark or gluon- 
jet with a final state formed by a cluster of strings corresponding to a given value of r]. 
We will follow the reference [65]. In Fig. 10 we show the formation of several strings in 
a nucleus-nucleus collision. At high density there will be overlapping of strings in the 
transverse space, forming clusters. Above the percolation critical density, essentially 
one cluster will be formed through the whole collision area. This colorless cluster can 
act as a non thermalized quark gluon plasma, where the color would no longer be 
confined to hadronic or flux tube dimensions but to the whole available area of the 
scattering. 

The string-like shape of the color fields are oriented along the collision axis. Con- 
sider a hard parton-parton collision, which produces two hardly scattered partons which 
are moving in the transverse plane to the reaction axis and in opposite directions (Fig. 
11). The interaction of these partons with QCD fields of the strings affects the parton 
momentum distribution, which determines the distribution of secondary particles in 
the jet. In particular, bremsstrahlung gluon radiation takes place in a tangential di- 
rection by the parton crossing the gluon field of each string. This gluon radiation will 
eventually produce low pt particles. In other words, the jets are not going to disappear 
but to be degraded to lower particles. This is at the origin of the difference found 
by STAR collaboration [66] between near and away side pt distributions. 

This interaction produces an asymmetry, the projection of momentum of secondary 
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particles in the jet along the reaction axis (Fig. 10, axis z) becomes larger than in 
the transverse direction (axis x). Let us considerer a quark qi passing through the flux 
strings transverse to the string axis. We assume constant static uniform chromoelectric 
field pointed along the reaction axis. Because in the static case chromomagnetic fields 
vanish, the force acting on the quark qi can be written as 

2 

where g in the QCD couphng, ^ — and A" are the Gell-Mann matrices. First, we 
compute the change of the transverse momentum of qi after crossing one string. Let 
axis X and y be the coordinate system in the section where crosses the string, let 
axis y be pointed along qis direction of motion and let ^ be the distance between axis 
y and line of qi movement. Then qi will cross the string surface in points {CuiiO) 
(C) 2/2(0) momentum change will be 



A pi (0 = / F (e, t)dt = F (e, y)dy = g^E'' (C, y)dy . (32) 
Doing the average over all ^ we obtain 

where Si is the transverse string area. Applying the Gauss theorem we have 

A p;= f ^^^^ f ^ A^A^ (34) 

^ 2ro 2 Jsi 2ro 2 Jv 2ro 8 ^ ^ ^ ^ 



where ^ is the color charge of quark q2 (string end) and z is the unit vector in the z 
direction. Doing the average over the color states of qi and q2, we obtain 

Ap, = ^^1^ = . (35) 

The average number of strings crossed by the quark is 

N = 2roLm (36) 
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where L is distance traveled by the quark, L ^ Ra and m is the number of strings per 
unit area 



hi A -\- A coUisions at mid rapidity, Ns ~ 2N'^^ at RHIC energies and SiV^^ at 
LHC energies. In the forward rapidity region Ng will be 2Na and 3A^^ at RHIC and 
LHC energies respectively. 

The mean transverse momentum of quark qi after it crosses N strings, would be 
< Ap|^ >=< {/^PZN-I + ^Pzf >=< > + < > +2 < /\pzN-l^PZ > ■ 

(38) 

As Ap^jv-i and /\pz may be pointed in any direction along the string independently, 
the last term vanishes and therefore 

< Ap >tot^ VNApi . (39) 

For as ~ 0.3 and Tq = 0.2 - 0.25, we obtain from (35) Api ~ 0.9 - 0.75 GeV/c. For 
central Au + Au coUisions at mid rapidity we obtain ^/N ~ 5.5 — 7. Then < Ap >tot— 
4.9 — 5.2 GeV/c, which is comparable with the pr triggered, therefore the back to 
back jet-like correlations disappears. On the contrary, for central d + Au collisions, 
Vn ~ 0.75 — 0.85 and < Ap >tot— 0.67 — 0.64 GeV/c, which is much smaller than 
the pt of the triggered jet, > 4 GeV/c, and the back to back jet-like correlation 
survives. An intermediate situation is Au + Au central collisions at forward rapidity, 
where ~ 2.5 - 3.2 and < Ap >tot- 2.2 - 2.4 GeV/c. In this case, the back to 
back jet structure is only partially destroyed. In the case of peripheral collisions, the 
suppression would be stronger when the jet is perpendicular to reaction plane (the plane 
spanned by the beam axis and the impact parameter, b ^ 0) than when the jet is in the 
reaction plane, in agreement with experimental data. Notice that before we deduced 
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stronger suppression at forward than at mid rapidities for the nuclear modification 
factor but we predict less suppression of the jet. In the first case, the origin of the 
further suppression is the normalization of the clustering in the initial state and in the 
second case, the lower interaction of the quark jet with the final state is due to the 
smaller string densities present at forward rapidities. Finally, let us mention that, in 
the framework of the CGC, it has been also found stronger depletion of the back to 
back correlations in p + ^4 and A + A collisions than in the case oip + p coUisions [67] . 

VI. CONCLUSIONS 

We have obtain a universal transverse momentum distribution that allows us to 
describe the low pr shape of the different particle species and the suppression of the 
high pt yield compared to the scaling with the number of collisions, Ncoii, expected on 
the basis of factorization theorem of QCD. 

The shape of this distribution is determined essentially by two functions, F{'r]) 
and k{r]), which depend on the density of color sources rj. The first one stems from 
the dynamics of the color clusters, and it is at the origin of the suppression of the 
multiplicity. It also controls the transverse momentum distribution of the produced 
particles due to the fragmentation of the cluster. The second one, k{r]), is related to 
the fiuctuations in the number of strings per cluster. At low density, k decreases with 
the density and, on the contrary, at high density it increases with the density. 

The fact that for p + Au and d + Au central coUisions the densities are not far 
from the minimum of k explains the Cronin effect. It will disappear at higher energy 
or higher densities because of the increase of fc. In the forward rapidity region, the 
number of color sources scales with the number of participant nucleons, Nai while 
in the central rapidity region it scales with the number of collisions Ncoii oc . 
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Therefore, comparing the pr yield scaled by the number of collisions, in the forward 
rapidity region there will be a suppression due to the additional factor 1/N\^^. 

The multiplicity distributions are related to the transverse momentum ones through 
the cluster size distribution, whose width is controlled by k. At high mean multiplicities 
the relative dispersion of the multiplicity distribution is given by 1/k, decreasing as 
the energy increases. 

The disappearance of the back to back jet-hke hadron correlations in Au + Au 
collisions is a final state interaction due to the interaction of the quark-jet with the 
cluster of strings formed in the initial state. This cluster is less dense in the forward 
rapidity region, so we predict smaller suppression of the correlation in this case. 

We can conclude that the clustering of color sources provides us a framework suit- 
able to describe the main features of the transverse momentum and multiplicity dis- 
tributions at RHIC and SPS energies. This, together with the successful explanation 
of the dependence of the fluctuations of transverse momentum on the number of par- 
ticipants and the dependence on the multiplicity of the strength of the two-body and 
three body correlations, favors this approach. 

We are aware of the limitations of our framework, coming mainly from the lack of a 
direct derivation from QCD. Also, we do not included any hard component in the single 
string. Therefore, at fixed energy and above some prm our picture will fail. However, as 
the energy and the density grow, there are more overlapping strings which extend the 
validity of the description to a higher pTm- Finally, asymptotically, our picture is valid 
for all transverse momentum, recovering a single exponential, exp{—F{r))p'^/ < p^ >j). 

Our approach has similarities with the CGC. In both approaches there is clustering 
-of gluons in the CGC and of color sources in our approach-. In both approaches, the 
initial state interactions -gluon saturation in the CGC or clustering of color sources 
in the percolation approach- produce a suppression of the px distributions. In CGC 
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the gluon field is renormalized as the color charge increases. In our approach, we have 
to redefined our main variables for the cluster as the color density increases. The 
suppression of the high pt yield and the saturation of the multiplicity per participant 
at high density are related to each other and they arc a consequence of the saturation 
of gluons at the scale given by the saturation momentum Qs in the CGC [68]. The 
clustering of strings produces also the suppression of pr and the independence of the 
multiplicity per participant on Na- We obtain, as in the CGC, larger suppression 
of the Pt yield at forward rapidity, and also we predict that the Cronin effect would 
disappear at higher energies. Both approaches obtain a scaling for the transverse 
momentum distributions. We do an extension of the soft spectra to the region of hard 
Pt by means of clustering of strings. The CGC does an extension of the hard region 
to the soft one by means of clustering of gluons, therefore in a broad range of pt both 
approaches should coincide. 

In QCD, the distribution for n gluons is approximately given, in the modified leading 
log approximation, by the generalized gamma function [69] 



whose width becomes broader as the energy increases. According to our considerations 
we expect a narrower width of this distribution above a certain density. It would be 
welcome the proof of that. 

Let us remark that other possibility of testing our approach is to study the behaviour 
of the forward-backward correlations which are proportional to the fiuctuations in the 
number of independent clusters. At moderate density, we expect a rise of these corre- 
lations as it is observed ai p + p and h + A collisions. On the contrary, at high density, 
these long range correlations should disappear as a consequence of the formation of a 
large cluster of strings [70]. Similar behaviour is expect in the CGC. 




(40) 



27 



Finally, let us comment that our approach provide us with an initial state which 
can affect the final state interactions, as it was shown above when we studied the 
disappearance of back jet-like correlations. 

This means that other phenomena as jet-quenching or interaction with comovers 
are not necessarily contradictory to our approach. We only claim that at low and 
intermediate pr, Pr < 4 — 5 GeV/c, the percolation of color sources is able to give a 
simple and reasonable description of data. 

We thank M. A. Braun, N. Armesto, F. del Moral, C. A. Salgado for dis- 
cussions. This work has been done under contracts FPA2002-01161 of CI- 
CyT of Spain, POCTI/36291/FIS/2000 of Portugal, Feder funds from EU and 
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FIGURES 




FIG. 1. Schematic representation of the number of clusters as a function of the number 

of strings of each cluster at three different centralitics (the continuous line corresponds to the 
most peripheral one and the pointed line to the most central one) 
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FIG. 2. Dependence of k on 77, using eqs. (23) and (24) (dotted line) and changing 
r) — 3.3?7 (solid line). 
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FIG. 3. Experimental PHOBOS data on low px distributions for pions, kaons and protons, 
together with our results. 
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FIG. 4. Comparison between our results and experimental data from Au-Au central and 
peripheral collisions and p + p collisions at s/s = 200 GeV. 
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FIG. 5. Comparison between our results and experimental data on the nuclear modifica- 
tion factor from Au-Au central and peripheral collisions at -s/s = 200 GeV. 
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FIG. 6. Our results for d-Au minimum bias collisions at -y/s = 200 and mid rapidity 
compared to experimental data. 
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FIG. 7. Our results for d-Au minimum bias collisions at -s/s = 200 and mid rapidity nor- 
malized to our results for p-p minimum bias collisions multiplied by the number of collisions, 
compared to experimental data. 
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FIG. 8. Nuclear modification factor for Ph + Ph central collisions at SPS energies. 
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FIG. 9. Predictions for Pb + Pb central collisions at 5500 GeV. 
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FIG. 10. Formation of strings in a nucleus-nucleus collision. 
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FIG. 11. Hard scattered partons that can interact with the strings. 
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